Molecular transport junctions: Current from electronic excitations in the leads 
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Using a model comprising a 2-level bridge connecting free electron reservoirs we show that coupling 
of a molecular bridge to electron-hole excitations in the leads can markedly effect the source-drain 
current through a molecular junction. In some cases, e.g. molecules that exhibit strong charge 
transfer transitions, the contribution from electron-hole excitations can exceed the Landauer elastic 
current and dominate the observed conduction. 
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Introduction. Electron transport in molecular tun- 
nel junctions has been the focus of intense recent re- 
searcbiiSi^. Theoretical modeling of tunnel conduc- 
tion^ starts from Hamiltonians that contain electron 
transfer (tunneling) interactions between molecule and 
leads as essential elements for current transport in such 
junctions. At the same time, energy transfer interac- 
tions — excitation/de-excitation of the molecule accom- 
panied by electron- hole (EH) pair annihilation/creation 
in the metal — are known to strongly affect the life- 
time of excited molecules near metal surfaces'". An essen- 
tial difference between these interactions is that electron 
transfer is a tunneling process that depends exponentially 
on the molecule-metal distance, while energy transfer is 
associated with dipolar coupling that scales like the in- 
verse cube of this distance, and can therefore dominate 
at larger distances. 

How will such dipolar interactions affect the conduc- 
tion properties of molecular junctions? Here we address 
this question by using the non-equilibrium Green func- 
tion (NEGF) formalism to derive an expression for the 
conduction in junction model that contains both elec- 
tron and energy transfer interactions, then analyze sev- 
eral examples with reasonable parameters. We conclude 
that current caused by electron-hole excitations in the 
leads may be significant, sometimes even dominant, in 
situations when strong asymmetry, of a particular type 
explained below, in the molecule- lead coupling is present. 
As a simple extreme example consider the case where the 
highest occupied molecular orbital (HOMO) is coupled 
only to one lead, while the lowest unoccupied molecu- 
lar orbital (LUMO) is coupled only to the other. Such a 
junction cannot pass current (in the absence of electronic 
correlations) if the dipolar interaction is absent. Realis- 
tic situations will not be that extreme, still whenever the 
HOMO-LUMO transition is of the charge transfer type, 
we expect some degree of such asymmetry. We have re- 
cently shown* that such situations may give rise to light 
induced current under zero voltage. We show below that 
also the current- voltage characteristic of such junctions is 
strongly affected by dipolar energy transfer interactions 
between molecule and leads. 

Model and Method. We consider a tunneling junction 
consisting of a molecule positioned between two metal 



contacts {L and R). The molecule is represented by its 
highest occupied molecular orbital (HOMO), |1 >, and 
lowest unoccupied molecular orbital (LUMO), |2 >, with 
energies £i and £2 and gap £21 = £2 — £i- The contacts 
are assumed to be free electron reservoirs, each at its own 
equilibrium, characterized by electronic chemical poten- 
tials /zl and hr, where the difference — /i/?, ~ is 
the imposed voltage. The corresponding Hamiltonian is 



H = Hq + Vm + Vn 
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where H.c. denotes Hcrmitian conjugate. Here the opera- 
tors Cm and cIj^ are annihilation and creation operators of 
electrons in the bridge (m = 1,2), while Ck and are an- 
nihilation and creation operators of electrons in the leads. 
The Hamiltonian Hq is a sum of terms that correspond to 
the isolated molecule (represented by its HOMO-LUMO 
levels in our model) and free contacts. Vm describes the 
electron transfer (tunneling) process between these sub- 
systems. This is the term usually employed to treat cur- 
rent in the biased junction. Vm represents coupling of 
the molecular HOMO-LUMO transition to electron-hole 
excitations in the contacts and is often used in models of 
energy transfer between the molecule and the contacts. 

In the Keldysh NEGF formalism,^ the steady-state cur- 
rent through the junction is given 

Isa = ±lj ^Tr [S< k(£;)G>(£;) - 11>,^{E)G<{E)] 

(5) 

calculated at the left {K = L with sign) or right 
[K = R with "— " sign) contact, where the direction from 
left to right chosen positive. 

The lesser and greater Green functions, G^^^, needed 
in jSJ can be obtained from the Keldysh equation 



G<'>(£:) = G""(£:) S<'>(£;) G'^iE) 



(6) 
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where the retarded and advanced Green functions, C"'", 
are given by the Dyson equation 

G'iE) = [E-H.^- ; G^(£;) = [G'"(£;)]^ 

(7) 

Here H™ is a matrix that corresponds to the molecular 
part (first term on the right) of the Hamiltonian (0) and 
is the retarded self-energy matrix due to both di- 
rect and dipolar coupling to the leads. Here and below 
the matrices are given in the basis {|1 >, |2 >}. 

The self-energies needed in Eqs.Q-JTI) are obtained 
within the usual diagrammatic technique on the Keldysh 
contour. In the non-crossing approximation (NCA)^^ 
this leads to 



(8) 



On the Keldysh contour these self energies areiSii^ 2x2 
matrices in the bridge space 
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Here and below we use S^^^^ to denote mrn' matrix 
element {m,m' = 1,2) of the self energy, P — M ot N 
corresponds to the physical process (electron or energy 
transfer to the metal, respectively) and K = L or R 
denotes the left and right leads, respectively, gk is the free 
electron Green function in state k, and fh = 2Sm^i+Sm,2, 
i.e. m = 1 if m = 2 and vice versa. 

After projection onto the real time axis we get the 
retarded, advanced, lesser, and greater components of 
these self-energies, which in steady state situations can 
be expressed in ener gy s pace. In the wide-band approxi- 
mation (see e.g. Ref.flir) the self-energies associated with 
electron exchange between molecule and leads have the 
familiar forms ^-S^ 
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(we neglect level mixing due to coupling to the contacts) 
where fiK is the chemical potential of the leads and K — 
L, R denotes the left and right electrode, respectively. 

The Langreth projection rulesi^ give the lesser and 
greater projections of the self-energies due to electron- 



hole excitations, Ea. (|10|l . in the (diagonal) form^ 

= I^B(^\cj,t,K)G<^{E + u;) (12a) 
SL^^)>(i?)= f^B^^\co,^K)G>^iE-co) (12b) 
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Here ''(i^) is density of electron-hole excitations in the 
lead K, p'^-^iLo) = SdEC^^\E,u:)fK{E)[l - fK{E + 
u;)], with C^^HE,cj) - Ek^k'eK^iE ~ ek)S{E + u; ~ 
ey) K, p[E)p{E + tj), where p{E) is the density of lead 
electronic states. In the spirit of the wide-band approxi- 
mation one can assume C'^^-' to be constant. This leads 
to p'K^iuj) = oj C(^) / [l - e-^'^] . Below we will use this 
expression in Ijl2dp to get B^^^ {lj, pk)- The retarded and 
advanced self-energies, ^ are difficult to calculate 

from the Langreth rules. For simplicity we assume, in the 
spirit of the wide band approximation, that all diagonal 
components of S^i^^) r ^^^^ ^.(nk) a(^) ^ ^^(nk) r(^)j t 
are purely imaginary. ThenS 

5,(NK)r(^) ^ [S(NK)>(j;)„s(NK)<(^)] ^ _ 

^ ^ (13) 

Eqs. ©,0, ((najl . (|12b|) and ^ have to be solved 
self-consistently until convergence is achieved. We use 
the level populations, = — * / ^G^^{E) (m = 1,2), 
as a test for convergence. Convergence is declared 
when the population values at subsequent iteration steps 
change no more than the predefined tolerance, taken be- 
low 10^6. 

The additive structure of the self-energy, Eq. ^ , makes 
it possible to separate the lesser and greater Green func- 
tions, Eq.©, and consequently also the source-drain cur- 
rent, Eq.Q, into contributions due to direct electron 
transfer to the leads and coupling to the electron-hole 
excitations. 



(14) 



is the usual Landauer expression for elastic current 

jL ^ ^ dE_ y{ML)f^r ( p\-p(M R) f-ia ( p\ 

sd f- / o_ / ^ m ^mm\ ) m ^ mm\ J 



m=\a 
X UdE) - fniE)] 



(15) 



while I^^ ' is the contribution from the electron-hole ex- 
citation. A simple expression for this current can be ob- 
tained when v'^m^^ <^ £21, where S^^^^*^'-^ takes an 
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approximate form 
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Here B^^^ is assumed to be a constant. Note that using 
B^-^") as a parameter in the full self-consistent calcula- 
tions below means that in fact we take |2(7(^) = 
B'^^'^ /2TTe2i- Using (|16() in jSJ leads to the electron-hole 
excitations part of the flux in the form 

(17) 



(ML) 



where 



^ rn I -L rn i ri 

-*/|f \Gl,jE)\'Y.'il^^<{E) (™ = 1,2). 

Further simplification of Eqs. (|15|l and (|17|l is achieved 
for strong bias, e.g. for negatively biased left electrode 
where /il ^ £2 and jin <C £1 so that fL = ^ and — 
in the energy range relevant to the integral in Ea. (|15|l . 
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(m = 1,2) cam be used in p7|) . Similar consideration 
apply in the opposite bias, leading finally to 



E 

m=l,2 



-p(ML)p(M_R) 



sgn(^L - A^fl) 



(18) 
(19) 



F 



(ML)-p(M_R) 
^ 1 



F1F2 



{ML)-^{MR) 



F1F2 



-^(^i^ - Ml) 



where is the step function d{x) = 1 for a; > and 
for a; < 0. Note that can be asymmetric to bias 

reversal (see also Fig. ^and discussion below). For, e.g., 
A*L > Mi? see that the magnitude of I^j'^ is deter- 
mined both by the value of i? = B'-^'' + B'-^^ and by 
the product p(-'^^)p(*'^^) t^ThHq /j^ ig determined by the 

product F^^^F^^^ (m = 1,2). The electron-hole con- 
tribution to the source-drain current is significant when 



^(ML) 



> F 



(MR) 



and/or F^ 



(MR) 



Below we com- 



pare the magnitude of the two contributions to the cur- 
rent for different junctions parameters. 
Numerical results. In the calculations reported below 
we used the following 'standard' choice of parameters: 
T = 300 K, £1 = eV, £2 = 2 eV, F^*^^ = F^*^^ = 0.2 eV. 
Values of other parameters are indicated in the figures. 
The Fermi energy is taken in the middle of the HOMO- 
LUMO gap. Chemical potentials in the left and right 
leads are assumed to shift with the voltage bias symmet- 
rically relative to the Fermi energy. Numerical integra- 
tion was done on the energy grid spanning range from 
-3 to 5 eV with the step lO'^ eV. 
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FIG. 1: The source-drain current Isd vs. applied voltage 
Shown are the total current I^d (full line, red) as well as 
contributions due to direct electron transfer jf^ (dashed line, 
green) and electron-hole excitations I^^'' (dotted line, blue) 
t(ml) _ p(Ji^fl) — Q I gY (a) ajjjj asymmetric 



for symmetric r\ . 

^(ML/R) ^ 

(b) cases. 



-.(ML) 



0.19 eV, and F 



(MR) 



0.01 eV 



Figure ^ depicts the current-voltage characteristic of 
the junction for the cases of symmetric and asymmet- 
ric coupling between the molecular LUMO and the con- 
tacts. Shown are the total current and its two compo- 
nents. In the symmetric case the current is dominated 
by the usual elastic electron (hole) transport through the 
LUMO (HOMO), and is symmetric relative to voltage 
reversal. The asymmetric case shows a significant contri- 
bution of the current associated with electron-hole exci- 
tations. The following points are noteworthy: (1) Z^^'' is 
significant when the LUMO is coupled asymmetrically to 
the two electrodes (2) This effect is particularly strong 
when the LUMO is coupled more strongly to the neg- 
atively biased electrode (i.e. <& < when the LUMO 
couples strongly to the left). Indeed, I^^^ is expected 
to be pronounced when the LUMO is populated and the 
HOMO is empty, which happens at such bias. Note also 
that the total current is asymmetric relative to bias re- 
versal in this case. 

Figure [3 shows the results of a model study 
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FIG. 2: Source-drain current Isd at $ = 3 V vs. molecule- 
contacts distance R. See text for choice of the coupling pa- 
rameters. Shown are the total Isd current (full line, red) 
as well as contributions due to direct electron transfer /f^ 
(dashed line, green) and electron- hole excitations /^j'' (dot- 
ted line, blue). 



of the dependence of the source-drain current on 
the molecule-lead distance R. We take Tm^^^ = 
Am^'^^ exp —a'm^^R to reflect a tunneling transition, 

while B'^^^R) is assumed to have a dipolar distance de- 
pendencei, i?^^^ = /3^-^^/R^. The parameters used are 



A. 



{ML) 



0.52 eV, A^2^^^ — nance of /^^ at large molecule-lead separations. 



and 



K) 



0.01 eV A3 



The choice of Ai^^^^ re- 



AY""-' = A'f'"^' = 0.27 eV, 
0.027 eV, ai!:^^^ = 1 
{K ^ L,R and m = 1,2) 
fleets a total lifetime broadening for electron transfer 
into the electrodes of 0.2 cV at a distance (from each 
electrode) of 1 A. The choice of /S*-^^ corresponds to 
taking B'^^^ = 0.01 eV at this distance. The relative 
importance of the and /^J''* components of Isd de- 
pends on the details of the molecule- leads couplings. In 



particular, when i? ^ oo for both leads Tm^^ /B 
(to = 1,2). In this limit I^^^ can become larger than 
/j^. As a specific example consider the situation where 

j,(M) ^ j.{M) ^^^^ j^^^^^ ^ p(ML)^-j,CA/)^ rj,^^^ ^gjj^g 

Eqs. ©, (Cni, (El), and for the case fj.R ei < 
£2 < ItL we obtain I^a'^/Isd ^ ?2/a(l " Ci - 6), oo 
and (1 - - 6)|1 - a - 61 when 6 < 1 ^ 6, 

^2 = 1 — Si and ^2 > 1 — Cij respectively. Our choice 
of parameters in Figure |21 corresponds to the third case 
and yields ultimate dominance of I^^'^ at large distances 
with I^J^/Isd ~ 100 as i? ^ oo. Note that this limit- 
ing behavior is obtained only when both left and right 
molecule-metal couplings decrease together. Experimen- 
tally one of these distances can be controlled by moving 
a tip while the other can be changed by adding insulating 
layers between molecule and substrat o^^'^7 . 
Conclusion. We have studied, within a simple model, 
the effect of dipolar energy-transfer interaction between 
molecule and leads on molecular conduction. We found 
that such interaction, that leads to electron-hole excita- 
tions in the contacts, can affect the current voltage char- 
acteristic of the junction in a substantial way and can 
not in general be disregarded. The contribution of this 
interaction can dominate the overall conduction for par- 
ticular asymmetric coupling where the molecular LUMO 
and/or HOMO are coupled differently to different leads. 
In addition, because of the different dependence of elec- 
tron and energy transfer on the molecule-leads distance, 
the relative importance of I^^ and I^^'^ depends on this 
distance, and can, in some cases, result in strong domi- 
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